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ABSTRACT 


The object of this study was to determine the feasi- 
bility of automatically controlling propeller pitch to 
reduce RPM fluctuations caused by variations in propeller 
speed of advance. The investigation was made ignoring the 
reduced frequency effect and power plant governing. 


A Victory Ship (AP2) was chosen as a prototype and a 
feedback system controlling the pitch of a controllable 
pitch propeller was developed. A non-linear analogue 
computer program approximating the system and the propeller 
Characteristics is presented. Electronic difficulties 
prevented obtaining total system response from this progran. 
The computer study indicated that a linear analysis would 
be useful. 


The linear analysis indicated that, for the prototype, 
in order to achieve a reduction in RPM fluctuation of 10% 
of that present in the fixed pitch case, a hydralic pump 
consuming over 300 HP would be needed. Systems utilizing 
RPM, shaft torque, and thrust are analyzed. Thrust sensing 
provided the poorest response and in fact increased the RPM 
fluctuation for lower frequencies. 


RPM governing by automatic control of propeller pitch 
was found to be inadvisable unless major design changes 


are made to reduce the spindle torque of presently available 
controllable pitch propellers. 


Thesis Supervisor: S. Curtis Powell 


Title: Associate Professor of Marine Engineering 
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I. 
INTRODUCTION 





The propulsion systems of ships operating in a seaway 
undergo speed fluctuations caused by variations in propeller 
load. The shaft torque produced by the propeller is a 
function of the flow velocity across the propeller disc 
area. As the ship rolls, pitches, and encounters waves 
this flow velocity fluctuates causing a fluctuation in 
shaft torque which the propulsion system must overcome. 

The propulsion system responds to the fluctuating load by 
allowing fluctuations in RPM. The amount of RPM fluctuation 
present depends on the system inertia, damping, and the 
effectiveness of the governing, if provided, of the power 
plant. High speed turbine installations with large 
reduction gears experience less RPM fluctuation difficulties 
than do direct drive diesels.” 

The largest variation of flow velocity averaged over 
the propeller disc area, henceforth called speed of advance, 


will probably occur at resonant pitching. Model tests 





* Evans, Harvey J. Conversations before seminar group 
of Department of Naval Architecture and Marine Engineer-~ 
ing, M.I.T. 24 April 1962. 
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indicate that at this condition the variation of speed of 
advance could be as high as 0.43 of the steady state 
value. For a geared turbine ship such as the Victory 
Ship (AP2) with a fixed pitch propeller a fluctuation of 
about 6 RPM will result at the frequency of resonant 
pitching, which is about 0.03 cycles per second. Oalcu- 
lations leading to these figures are given in Appendix E. 
For systems with less effective inertia and damping the 
fluctuations would be correspondingly higher. 

Ships propellers are designed to operate efficiently 
at a designed speed of advance. Their efficiency is 
reduced when subjected to a constantly varying speed of 
advance. The propulsion system could deliver full power 
when the ship was operating in a seaway if the load could 
be made sufficiently constant. 

A controllable pitch propeller has the ability to 
alter the propeller pitch at any time. If the pitch were 
altered to counteract the effect of the changing speed of 
advance, a constant torque load would be provided to the 
power plant. The efficiency of the power plant and 


propeller might be increased and the entire system could 





## Burtis, Evenson M. and J.D. O'Connell. Measurement 
Of Ship Model Wake Variation in Rough Water, Naval 
Engineers Thesis, Department of Naval Architecture and 
Marine Ingineering, M.1I.T. Cambridge, 1959. 
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be operated in heavy seas at high power. A literature 
survey has revealed that no reported investigation of 
this possibility has been made. Extensive use of one 
governor systems controlling RPM by altering both pitch 
and throttle setting has been made by the aircraft 
industry in conjunction with turboprop aircraft. 

The purpose of this paper was to investigate the 
possibility of governing propulsion system RPM by auto- 
matic control of propeller pitch. Systems sensing RPM, 
torque, and propeller thrust were investigated. System 
speed governing by controlling the load may or may not 
be more advantageous than conventional governing of the 
power plant. The regulation of some variable other than 
speed might also provide a more optimum overall operation. 
This investigation, however, examines the feasibility of 
speed governing through propeller load control only. 
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Il. 
GENERAL PROCEDURE 





Several major analytical assumptions were made at 
the outset. Unsteady flow was trated using the quasi- 
steady approximation and the reduced frequency effects 
were neglected. The speed of advance was assumed to be 
a steady state value with a sinusoidal fluctuation super- 
imposed thereon. The control systems investigated 
contained a simple one loop feedback with no special 
compensating networks. The power plant was assumed to 
deliver constant torque and no model of the governor was 
developed. 

The standard method of presenting propeller character- 
istics envolves a plot of non-dimensional thrust and torque 
coefficients as functions of advance ratio, J, and pitch 
ratio me These plots are very nearly families of straight 
lines. For the non-linear study the propeller curves 
were assumed to be families of straight lines. The slopes 
and values of the actual and the assumed curves are the 
same at the operating point. The accuracy of the approxi- 
mation decreases as conditions depart from the operating 


point. The approximations lead to expressions (6) and (13) 





for propeller torque and thrust respectively as derived 


in Appendix A. 
Qp = Kypw* - KyoVyoW - Kz we eeee(6) 
T = Kap ww = K5VQW ~ Ke Ww ° (13) 


The K's are constants for a given propeller and 
depend on both the propeller curves and the dimensions 
of-the propeller. 

For the linear analysis, the propeller torque and 
thrust were assumed to be linear functions of speed of 
advance, pitch, and RPM. This is the customary incre- 
mental linear analysis and is valid only for small 
variations of the variables about the operating point. 
The equations for incremental variations in propeller 


torque and thrust are derived in Appendix E and are 


O 
hj@ 
i 


¥ 
Ky Vo = Ke + K*w (81) 
and 


& t 
KV, +K + K*w (100) 


3 
i 


The K's are found by assuming a small change in 
the corresponding variable (* 0.1), and while keeping 
all other independent variables constant, determining 
from the propeller curves the corresponding change in 
the dependent variable. fThe K™ then becomes the 
dependent variable change divided by 0.2 “the Troost 


Standard Series for Single Screw Merchant Vessels (B4-55) 
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propeller curves were used throughout. ™ 

No test data providing blade spindle torque directly 
ig available. The torque about the blade spindle axis 
was assumed to be composed of hydrodynamic, sliding 
friction, viscous friction, inertia, and centrifugal 
force components. As with the propeller torque and 
thrust coefficients, the hydrodynamic component of blade 
torque can be represented by a family of blade torque 
coefficient curves, Kp, where advance ratio and pitch 
ratio are the independent variables. These curves were 
also assumed to be families of straight lines and the 
resulting expression for hydrodynamic blade torque as 
derived in Appendix A is 

Qpy = K7pw® = KgVW - Kow? . (21) 

Reference (11) contains data onthe control rod 
forces for a controllable pitch propeller installation 
on a tug,(YTB-502), while underway at various speeds. 
This data also includes details of the pitch setting 
mechanism, ship's speed, centrifugal force components as 
functions of RPM and pitch, and other useful information. 
An approximate wake fraction of 0.164 was assumed based 


on the block coefficient of the tug. Since the tests were 





* Lewis, F.M. "Propeller Coefficients and the Powering 
of Ships," Society of Naval Architects and Marine 
Engineers-.Paper, Nov. 1951. 
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made under steady conditions, the only forces assume 
acting were hydrodynamic and centrifugal. After 
subtracting the centrifugal component, the control 
rod force was converted to a spindle torque, non- 


dimensionalized, and finally plotted as blade torque 


d 


coefficient versus advance ratio and pitch ratio. The 


results are shown in Figure B=-3 below and the calcu- 


lations leading to this plot are given in Appendix B. 


Torque tending to rotate blade t 
ahead pitch 
0 / ‘2 3 - Ss A 
= Vo 
nb 
FIGURE B=3 


Derived Kp Curves for YTB-502 
-7- 








As before stated, the Ky curves were assumed linear. 

The hydrodynamic spindle torque is a strong function of 
the location of the spindle axis relative to the hydro- 
dynamic center of pressure. It was assumed that the 
effect of moving the spindle axis would be to shift the 
Ky values up or down the family of Ky curves. Once the 
operating point of the prototype assumed for this study 
was chosen, the operating point of the Ks curves given 
in Figure Be3 was assumed moved to the prototypes 


simply by shifting the K, values corresponding to the 


B 
tugs operating point to the approximate position of the 
prototype operating point. This can be justified only 
in that such a propeller could be made. The slope and 
spacing of the Kp curves was assumed to remain the same. 


The centrifugal force component of blade torque as 


derived by Boswell and others ig” 


tip tip 
=. fos (lax ._ i in) 
2Be = PW N Axydr + 5 = Sin2d dr 
Hub root 


(28) 


This expression would be difficult to compute on an 


* Robert J. Boswell, A Method of Calculating the Spindle 


Torgue of a Controllable-Pitch Propeller at Design 
Conditions. David Taylor Model Basis Report No. 1529, 


Aug. 1961, pp. 31-34. 





analogue computer and a simplification given in Appendix 
A leads to expression (35) where the constants Kio and Ki, 
are functions of propeller geometry alone. 

Q@pc = = Kip Me = wep - (35) 

The sliding friction component of spindle torque was 
obtained from altering the usual expression for sliding 
friction in a journal bearing to allow for a moment 
tending to rotate the spindle shaft about a line perpen- 
dicular to the shaft centerline. The resulting expression 


is 
ert, AN 
QBs = sd = (36) 


The double sign in expression (38) indicates that the 
direction of the torque must be taken either (+) or (-) 
depending on whether the direction of pitch change is 
toward lower or higher ahead pitch respectively. 

The remaining components of blade spindle torque are 
obtained by direct application of their respective 


definitions. The final expression for blade torque is 
hydrodynamic 
Qa = K7pw)* - K8VoW - Kg w ® 


sliding friction viscous friction 
et. 
+ ert, f N+ LpW3N 
b 


inertia centrifugal 


OD > 


i. - 2 
+ JjpW ZN Kip OW - Ky wp (38) 
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There are numerous variations of hub mechanisms in 
use. All are required to convert motion in the direction 
of the propeller shaft axis into rotation of the blades 
about their spindle axes. As shown in Appendix A, each 
of the four mechanisms most frequently used result in 
approximately the same amount of blade rotation for the 
same input of linear motion along the propeller axis up 
to a pitch ratio of 1.5. The same is true for the force- 
to-torque transmission. For this reason the simplest 
mathematical relation of the four was used throughout. 
The results are approximately applicable to all the 


mechanisms. The expressions used are as follows: 


ae: 
w == (39) 
and 
1 


There are two general types of controllable pitch 
propeller mechanisms. The first, designated herein as 
Design I, utilizes a control rod extending through the 
propeller shaft and is driven by a hydralic piston with- 
in the ship. The second, designated herein as Design II, 
is driven by a power piston housed within the propeller 
hub. The power piston of Design II receives oil from 
&@ pump within the ship through oil passages in the 
propeller shaft. Design II, as presently available, 
includes a safety spring within the hub that makes the 
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propeller "Fail-safe" by forcing the blades into ahead 
full pitch in the event of a failure in hydralic oil 
pressure. 

The equation of motion of both these designs is 
obtained by equating the oil flow from the pump with 
the sum of the flows into eens et sn which include 
flow from piston motion, flow from leakage, and flow 
due to compression of the fluid and expansion of piping. 
The final equation of both systems are of the same form 
With the exception of additional terms in the Design II 
expression that are generated by a safety spring force. 
Design I was used for all further study. Schematics of 
both designs are given as Figures A-7 and A-8 in Appendix 
A. . The equation of motion of Design I as derived in 
Appendix A is 


Input Ram Flow Leakage flow 
Cas 
Kpx = iéRAp + Ky Kop W Kav Ww 
D AR R R 
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The Victory Ship (AP2) was chosen as a prototype to 
determine values of the various parameters. This is a 6000 HP 
geared turbine merchant ship for which considerable data was 
available. The detailed determination of the various 
parameters present in equation (42) above is carried out in 
Appendix B. The propeller blade was assumed to be an 
elliptical flat plate for the determination of blade inertia 
and centrifugal force components of blade spindle torque. 
The dimensions of the hub were obtained by scaling up the 
design given by Rupp in Reference ll. 

The three remaining parameters to fix were area of the 
ram, the size of the pump, and the hydralic pressure. Since 
the initial cost and maintenance on systems using more than 
1000 psi would be high, this value was chosen as the 
maximum hydralic pressure. The maximum diameter of power 
piston that could be accommodated within the propeller hub 
has a radius of 15.75 inches. 

Although Design I has the power piston within the ship, 
this was assumed the largest piston size that would be fitted. 
Calculations given in Appendix B show that with the assumed 


hydralic pressure of 1000 psi and this ram area there is still 


=e 





a likelihood of stalling the mechanism (by opening relief 
valves) even with this large size power piston. The 
maximum combinations of blade centrifugal, hydrodynamic, 


and sliding friction when acting together require 1000 in® 


of piston area. The assumed piston area is 775 in, 
Larger piston sizes would require larger pump capacities 
which are excessive for this piston area. 
As shown in Appendix E, the time constant of the pitch 
control mechanism is 
a = (60) 


Kp DM 


The effectiveness of the system in regulating RPM 
fluctuations requires that this time constant be small. 

The larger pump capacities correspond to larger Kp: Assuming 
that the marine industry would not accept a device consuming 
more than 5% of the installed main plant horsepower, a 
maximum pump size of 300 HP was assumed. Kp was then 
determined from manufacturers data. 

The viscous friction forces were obtained from an 
analysis of transient spin test data taken by Rupp. This 
data consists of a plot of control rod forces versus pitch 
for the operating conditions of going from full ahead to 
full astern and vice-versa. This data is reproduced as 
Figure B-7 below. The thrust forces causing high values 


of sliding friction were assumed zero for this test. The 
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YTB-502 Control Rod Force for Spin Test 
in Air From Data by Rupp 


long region of both transient curves that are parallel 

to the steady state condition are assumed to correspond 

to the condition of maximum pump output and constant piston 
speed. The deviation from the steady state curve in both 
cases is approximately 3 x 10 lbs. This control rod 

force was converted to spindle torque and scaled up by the 
ratio of propeller diameter squared since viscous friction 
ls proportional to surface area. The speed of rotation 
was assumed constant and equal to the total angle divided 


by total time. 
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The substitution of the various parameters into 
equation (42) revealed that several of the terms would 
be very small in comparison to others. The coefficients 
of each independent variable was examined and those 
terms having coefficients of less than 1% of those re- 
maining for the same variable were eliminated. A 
further simplification was made by setting to zero all 
terms including the first derivative of shaft speed. 
This assumption is based on the knowledge that propeller 
RPM does not vary to a large extent even without regu- 
lation and should vary less with regulation. The 
resulting simplification of equation (46) is 


. K7-K e 
Kpx = p Bee. ms Al pi ~ ao teeS 
WD Ap R R 
~ (Kg#Kiglw? + (fp + fri f N 
R bR 
K ° 

+ sls - KWo +) (oC BY + Wp) (65) 
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EEE 
PROCEDURE AND RESULTS OF 
THE NON-LINEAR ANALOGUE COMPUTER STUDY 


A. Procedure 

Equations (6), (13), and (65) given above envolve 
non-linearities in the form of products of independent 
variables and a discontinuous sliding friction term. [If 
saturation of the pump is also allowed, another discon- 
tinuity is introduced which would limit the value of pump 
output, Kpx, to some finite value. All of these non- 
linearities can be programed on the analogue computer 
that was available. The program includes four multipliers 
to generate ie pw, Vand Vow. It also includes a 
switching mechanism controlled by the derivative of pitch 
to approximate sliding friction. 

The program was arranged so that shaft torque and 
thrust were generated in order to test the final system 
sensing these variables. 

The program is shown in Figure C-2 below. The para- 
meters determined in Appendix B were then inserted into 
the program and voltage scales were assumed for each 


variable. The gains of various loops were then adjusted 
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to provide signals that would not saturate the various 
amplifiers and multipliers of the computer. The details 
of this scaling is given in Appendix C. 

The input to the system was taken to be the voltage 
analogue corresponding to speed of advance, Vo» and 
consisted of a steady positive value with a sinusoidal 
voltage of amplitude AV, provided by a signal generator. 
The output consisted of a voltage analogue for shaft 
speed, W , and both input and output were recorded on a 
two channel Sanborn recorder. It was hoped that data over 
a range of frequencies could be obtained in this manner 
for the fixed pitch case as well as the system operating 


in closed loop sensing shaft speed, torque, and thrust. 


B. Results of the Non-Linear Study 
The frequency response of the fixed pitch system was 


Obtained. The data is recorded in Appendix C and is 
shown graphically as a plot of the magnitude of the ratio 
of output to input in decilogs in Figure I below. The 
strong stabilizing effect of sliding friction was very 
apparent. When the signal corresponding to sliding 
friction feedback was removed with a sinusoidal AV, 
input of amplitude 6 ft/sec the second derivative of 
pitch, ~, increased by a factor of 80. The feed forward 
path containing Vo had no effect on p even for very 


large values of AV). 
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Values generated by the computer for propeller 
torque, thrust, and blade torque were within 5% of the 
calculated values. The time constant of the system with 
fixed pitch was observed by injecting a step in engine 
torque and recording the decay in shaft speed variation. 
The time constant was observed to be approximately 3.28 
seconds. The decay to steady state required about 10 
seconds (approximately 3 time ae 

Electronic difficulties, the details of which are 
given in Appendix D, prohibited obtaining data for the 
control system operating in closed loop. 

In all cases the node containing } was very small 


and consisted primarily of computer generated noise. 


-20- 





IV. 
PROCEDURE AND RESULTS OF THE LINEAR STUDY 
A. Procedure 
The results of the non-linear study indicated that 
valid further simplifications could be made in equation 
(46). Eliminating the terms containing Vw and >} as 


suggested by Chapter III above, equation (46) reduces to 


; L K 2 
P = RAR Kpx = au (K7-Kj1) ll (67) 
TD i AR R 


- KeVqw - (KotK10) (wo : (tpetrt sf N )| 
a “R bR 


Equation (67) includes ram flow and leakage terms 
only. Since the linear analysis is valid only for small 
variations around the operating point, the leakage was 
assumed constant and is equivalent to an effective 
reduction in pump output Kp/*. As shown in Appendix E, 
equation (67) eventually reduces to a formula whose 


frequency domain transfer function is 


p(s) — G/M 


y'(s) TpS+1 


where y'(s) = Laplace transform of the incremental 





mechanical input to the pitch control mechanism. 

Using the linear approximations for propeller 
response before mentioned, and assuming the power plant- 
propeller system to be a first order system containing 
inertia and damping, the complete systems for sensing 
RPM, torque, and thrust were drawn in block diagram 
form as shown in Figure II. 

The gains G,, and Gg for the shaft speed and shaft 
torque sensing systems respectively were chosen so that 
a damping ratio of 0.4 was assumed by the closed loop 
system. Higher values of loop gain reduce the low 
frequency response but cause higher peaking or ampli- 
fication at a higher frequency. This peaking is apparent 
at about 2 radians per second as shown in Figure E-6. The 
gain for thrust, Gms» was taken to be a very low value 
only to demonstrate that any positive value of Gm is not 


beneficial in regulating RPM at the lower frequencies. 


Be. Results 

System gain is defined as the ratio of shaft speed 
variation for a given variation in speed of advance, a. 
The system gain for each of the three linear systems 
pictured in Figure II was obtained and plotted by the 
standard methods. 

The results are shown in Figure E-6 where system gain 


is plotted in decilogs. 
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(b) Regulation by Shaft Torque Sensing 
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Linear Approximation to Complete Systems 
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DISCUSSION OF RESULTS 





As shown in Figure E-6, the fixed pitch response 
determined on the analogue computer was in good agreement 
with the linear analysis for frequencies below 0.10 rad/sec. 
It is of interest that the input levels used to obtain this 
response curve are extremely high. The minimum value for 
the fluctuating component of the speed of advance is three 
volts which corresponds to thirty feet per second. The 
steady state value of speed of advance is only 16.7 ft/sec. 
The maximum value of input reached 40 volts corresponding 
to 400 ft/sec. These high values of input were necessary 
to obtain an output above the noise level. The results 
therefore indicate that the linear analysis should be 
reasonably valid even for moderate fluctuations around the 
operating point. 

Had it been known that the terms containing ‘p and V in 
equation (65) would be insignificant, an analogue eliminating 
these paths would have been developed. As explained in 
Appendix D, hum in the p path made the obtaining of data 
for the closed loop system impossible. The simpler analogue 
would not have this low signal level path and might have been 





more successful. Unfortunately time was not available 
to reprogran,. 
The simplifications leading to equation (65) as 


shown in Appendix C, call for ignoring the term 


dtp + fr AN ) 
at bR 
This friction term represents an oil flow and since 
it is discontinuous at p = 0, corresponds to an impulse or 
hydralic shock. This hydralic shock would very likely be 
damaging to relief valves and other parts of the systen. 
The introduction of an accumulator to cushion this shock 
would introduce additional compliance into the system, but 
as shown in Appendices C and D compliance is relatively 
unimportant. There are no compliance terms remaining in 
the final linear analysis. 
The simplifications of expression (46) given in 
Appendix C show that inertia loading is relatively unimpor- 
tant. Design I has a larger inertia loading because of the 
relatively larger mass of the control rod compared to the 
smaller power piston of Design II. The area of the high 
pressure port delivering full flow to the power piston of 
Design II would have to be 9.1 square inches to match the 
output of the variable displacement pump at the same pressure. 
Since inertia loading is relatively unimportant and since 


providing a large oil port area would be a difficult engineer- 
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ing problem, there is little advantage to Design II for 
this application. 

The transfer functions used in comparing the various 
systems are stable for all frequencies. These transfer 
functions are based on numerous assumptions that neglect 
many of the lags and storage elements actually present. 
They also do not consider nonlinearities such as sliding 
friction, stKiction, saturation, and the nonlinearities 
present in the propeller characteristics themselves. 
There is therefore a possibility that the actual hardware 
might exhibit matural frequencies corresponding to 
resonances in the remainder of the propulsion systen. 

The lowest frequency of torsional vibration resonance for 
the Victory Ship (AP2), as measured and reported in 
Reference 16, is 28 cycles per second or 175 radians per 
second. 

With the restrictions of a maximum pressure of 1000 
psi. hydralic pressure, a maximum pump power of 300 HP, 
and the approximate dimensions of present day controllable 
pitch propellers, the improvement in shaft RPM regulation 
over a fixed pitch propeller installation using automatic 
pitch control corresponds to a reduction in RPM fluctuations 
by a factor of nine. The improvement in propeller and 
power plant efficiency would not be great enough for this 
small an increase in RPM regulation to warrant the 


installation of such a large and complex systen. 


=a 





Several improvements in the system analyzed might 
reduce the size of the equipment necessary or improve 
the system response. The raising of the hydralic 
pressure would reduce the size of the piston necessary 
and thereby reduce the capacity of flow mecessary from 
the pump. This would improve the system by reducing 
the important pump-piston time lag. Increasing the 
horsepower of the pump also improves the system by 
decreasing the pump-piston time lag. Compensation in the 
form of lead or lag networks in the feedback loop would 
delay the peaking present in the proposed systems and 
allow a higher closed loop gain and hence a higher degree 
of regulation. 

It is unlikely that the above mentioned changes 
would result in improvementsfin system performance to any 
large degree. If spindle torque could be greatly reduced, 
the size of the power piston could be reduced and the 
entire system could be designed with better response 
and be made smaller. As shown by equation (38), spindle 
torque consists of many components, each of them large. 
By increasing hub size the spindle bearing could be made 
larger and the lever arm of the hub mechanism increased, 
The application of non-friction bearings might reduce 
this factor although roller and ball bearings are presently 
considered unacceptable for naval service due to shock 


susceptability. If the propeller were designed to operate 
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at one RPM, the hydrodynamic and centrifugal components 
might be reduced by careful design. These significant 
changes might make automatic control of pitch to remove 
fluctuation in RPM feasible. 

As demonstrated by Figure E-6 there is no advantage 
in sensing either torque or thrust in lieu of RPM when 
the variable being regulated is RPM. The improvement 
in regulation offered by RPM sensing is 9.7 decilogs 
which corresponds to a reduction in fluctuation by a 


factor of 9. 
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Vel 
CONCLUSIONS 





For geared turbine or geared diesel propulsion plants 
used by the U.S. maritime industry, automatic control of 
propeller pitch to reduce wave excited RPM fluctuations 
would be excessively expensive in horsepower, space, and 
complexity for a small gain in RPM regulation, 

Significant reductions in dynamic blade spindle torque 
through design changes would be necessary to make this 


scheme practical. 





Vik. 
RECOMMENDATIONS 





Automatic pitch control investigation requires a 
knowledge of spindle torque. Very little is known in 
this area. It is recommended that future effort be con- 
centrated on the determination of dynamic spindle torque. 

A more useful approach to the problem of automatic 
pitch control would include the presence of engine 
governing with a view toward developing a system that, 
with one governor would receive a bridge signal and 
control pitch and engine throttle setting to achieve the 


thrust or ships speed desired by the operator. 





VIII. 
APPENDIX 





APPENDIX A 


Derivation of System Equations 


1. The Power Plant-Propeller System: 
Figure A=-1 represents a schematic of the propeller- 
power plant system. 


The equation of motion of this system is 


$ t 
Qe - Qp - Lpwe (Jip + Jip) & ; (1) 
where 
QE = developed engine torque, lb-ft 
Qp = absorbed propeller torque, lb-ft 
' 
Jip = effective polar mass moment of inertia of the 2 
propeller (includes entrained water), lb-ft-sec 
Uy 
Jip = effective polar mass moment of inertia of the 
engine and gearing referred to the propeller shaft, 
lb-ft-sec® 
cL = mechanical damping in the power plant and gearing, 
Qs = shaft torque, ftelbs. 
WwW = shaft rotational speed, radians / second 
UO = shaft rotational acceleration, rad/sec* ‘ 


The torque delivered by the engine is 

QE -L WwW - 
The shaft torque expression is 

Qg = Qp -Lyw- Ji ph : (2) 
Curves of Ko for most propellers are nearly linear over a 


wide range of adrance coefficient and pitch ratio. Figure A-2 


=a 








FIGURE A-l 


The Propeller-Power Plant System 
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FIGURE 2 


Assumed Linear Torque Coefficient (Ko) Curves 
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represents an approximate general set of Kg curves assuming 
that they are completely linear. 


The equations of these lines are 


fo 
5 = Mad + boKy + hy (3) 


Kg * Guz DS 


propeller pitch, ft. 


D propeller diameter, ft. 


mig = slope of the Kg curve at the operating point of the 
true Ky curve and assumed constant for this analysis, 
non-dimensional 


J = propeller advance coefficient = = 

Vo = speed of advance for the propeller, ft/sec 

n = propeller speed, rev./sec 

bg = ratio —— » non-dimensional, assumed constant 
ho = ”/D intercept of the K=O line. 

e = water mass density 1.024! 


When MQ» OQ» and Qo correspond to those at an operating 
point taken from actual propeller curves, the response of the 
propeller to changes around the operating point will be 
accurate within the limitations of a steady flow analysis. As 
the propeller operates farther from this point the assumption 


of linearity becomes less valid, but is considered sufficiently 


1. This unit for mass density for water was taken to be 
consistant with Troost's propeller curves used to 
determine coefficients of equation (3). 
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accurate for this study. 


Solving equation (3) for Qp, 


Qp = £ 92p4 ~ MgVo pnpt - nae np? (4) 
bQ La Jal 
Since n= 2, (5) 
Tt 


equation (4) can be written in the form 


Qn = Kypw® - KevoW- Kzw* (6) 
where 
A 
K, = OD (7) 
FZ tg : 
4 
R=", (8) 


2 TT bg 


and K, = ho PD? (9) 
3 a e 


Substitution of equation (6) into equation (1) completes 
the derivation of the propeller-power plant system equation. 


Qn - Ky pw + Kpvow 2 pW = (Jip + Jip) & (10) 
E7 1Ppw + KoVo + KzW™ -ARW = D+ E 


The thrust coefficient for most propellers can also be 
considered linear with respect to ?/D and J. Again this is 
accurate for small changes around an operating point and less 
accurate for large departures from this point. Figure A=-3 
represents the linear approximation to the actual K- curves. 

The procedure for obtaining thrust as a function of Vo,w, 
and p follows in the same manner as used in obtaining the 
slmilar expression for Qp. 

From the linearized K+ curves 


Pr = mJ +b, K. + he © (1) 
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Linear Approximation to Kn Curves 
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Linearized Hydrodynamic Blade Torque Coefficient Curves 
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Using the definition of K_, 
ln 





<T= Spee > = 
and solving for iia ° 
Ty = Kypw “ K.VoW - Ke oe (iS ) 
where 
K ¢D° 14) 
+" aT” . 
Ko = ——1e) , (15) 
21 Ans 
ae * 
amd Ke = ° 16 
6 ant he (16) 
With b = 2 , ee), 
TT 


3 
= 


= slope of K+ curve, 


and h.= 2# intercept of K, = O line. 
iF D Q 


2. Spindle Torque 

The equations of motion of any servo system are a function 
of system loading. Loading in this case consists of blade 
spindle torque. Torque about the spindle axis can be viewed 
as composed of hydrodynamic, sliding friction, viscous friction, 
inertia, and centrifugal force components. 

These components are derived separately. The equation 
for spindle torque is 

Qp = QpH + Qps + Qpy + Qpz + Que (18) 
where 


Qp = total spindle torque that must be overcome by the 
pitch changing mechanism, lb-ft 


hydrodynamic spindle torque, lb-ft 


QBH 


sliding friction torque, lb-ft 
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Qpy = viscous friction torque, lb-ft 
QBJ = Jip Wp, lb-ft 


Qpo = centrifugal force torque from the mass of the 
propeller metal, lb-ft. 


Very few curves of spindle torque are available. Hydro-~- 
dynamic spindle torque is a strong function of the spindle 
axis position relative to the center of hydrodymic 
pressure. Tests made at M.I.T. show that curves of Kp where 
Kp is defined as 


K,= BH (19) 
B 6 n2D> ‘ 


are very nearly linear functions of advance coefficient and 
pitch ratio just as is c@ As with Qy linear curves for K, 
are assumed as shown in Figure A-4, and the expression for 
Qpy is derived based on this assumption. 

For these curves a spindle torque that tends to increase 
pitch is considered a positive torque. 

Proceeding as before, 


+ h, (20) 


v= Mpd + dg Kp, 


and 


2 


Qpy = Ky pw = KgVo - Kyw (21) 


where 


Kv = D (22) 
1* _hep , 





Conversation with S. Curtis Powell, March 2, 1962, 
M.I.T. 
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ke= ‘Tae (23) 


and Kg = - C D ; (24) 


TT bp 


The spindle is not a journal loaded in the normal 
Manner. Propeller thrust acts at a considerable moment 
arm tending to rotate the spinde about an axis perpen- 
dicular to its centerline as shown in Figure A-5. 

if the bearing were normally loaded the friction torque 
would be 

fPr , lb-ft 
where 

f = coefficient of journal friction 

P = bearing load, 1b 


spindle radius, ft. 


= KH 
Nl 


total thrust on blade perpendicular to the spindle 
axis, lb. 


An approximate expression for spindle sliding friction is 


where 


yi = lever arm from center of thrust on the blade to the 
center of the spindle bearing, ft. 


b = spindle bearing length, ft. 

N = number of blades 

Sliding and static friction are assumed to be the same. 
Blade viscous torque can be expressed as 


QBV — ae Bw BN (26) 
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FIGURE A=-5 
Schematic Indicating Blade Sliding Friction ( 


Qna) 





FIGURE A-6 


The Relationship Between Pitch and Blade Angle 
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on = coefficient of viscous friction, lb-ft-sec. 
CO, = spindle rotational velocity, radians/sec. 


Blade inertia torque is 
= 1 ye 
QB = Jip Lp N ) (27) 


where 


Cy 
J 
t= 
iH 


polar mass moment of inertia of the blade about 
the spindle axis, lb-ft-sec¢@ 


spindle rotational acceleration, rad/sec* ° 


€ 
bo 
1 


The expression for blade spindle torque as derived by 


Boswell and others is 3 


ti tip 
Qac = -Parn \ A nyde ate = = T minysin 2a 28) 


2. 
hub root 
where 
e = mass dgnsify of the blade material, 
P  lbesec“/ft 

UY) = propeller angular velocity in rad/sec. 
A = blade section area, £te 
X,y = coordinates of the blade section centroid 


measured in the plane of the section from 
the blade spindle axis, ft. 


( radius of the section, ft. 


i = the minimum, moment of inertia of the blade 
MIN 4 
section, ft 





3. Robert J. Boswell, A Method of Calculating the Spindle 


Torque of a Controllable-Pitch Propeller at Design 
Conditions. David Taylor Model Basin Report No. 1529, 


Aug. 1961, pp. 31-34. 
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TMax = the maximum moment of inertia of the blade 
section, ft 


cA = pitch angle of the section from its major 
axis, and is a function of both r and p, 


ORC = centrifugal component of blade torque, positive 
tending to rotate pitch ahead., lb-ft 


Equation (28) is not readily adapted to analogue 
computation. Several simplifying assumptions are necessary 
to make it so. 

Let 2 be the pitch angle in radians at the seven 
tenths radius. This angle corresponds to the "pitch" of the 
blade as shown in Figure A-6. | 

A in equation (21) will be greater than =) near the 
root and less than O, at the tip. 


Assume C= O This assumes a constant pitch angle 


B° 
over the blade length. Equation (28) can then be written: 
tip tip J 
- QBo = pain aay dr + pusksin26, (umm Lm) 4" (29) 
P P é 
hub root 


From Figure A-6 is can be seen that 
ate = tan y) (30) 
1D B 


and in fe = tanut 


40 
B =D ) (31) 


whose expansion is 


3 
5 _Ff L (€ Sr (32) 


The expansion for Sin 2 oF is 
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Sin = 20. — 4) Eee oe (35) 


au 
Equating the first terms of these expansions, 
2p * Sin 20, (34) 
TD 


This assumption will introduce an error for this 
component of spindle torque that is sometimes conservative 
and sometimes not conservative since the total spindle torque 
is a sum of components with changing signs leading to a sum 
that also changes sign. The maximum error will occur for 


large values of a ° For the case studies, (ft) = ie 


For this case Ms 


2F = 0./ 05, 
Tr D 


and Sin 2 o= 0.644 , 


With these assumptions equation (29) reduces to 


Qn = “K, uF » Ko (35) 
where va 
p 
K = 6 
10 N ~ a (36) 
hub 
d 
\ tip 
K,, = TD = 
mt te | (Te - sae (37) 
hub 


Substituting equations (21), (25), (26), (27), and (36) 
into equation (18) ) 
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Q., = Ky pw” - KgVow - Kyu? + £ sun + <6 Nie 
Jib dy = KyQW - Ky wep. (38) 


Where N = number of blades 

The double sign before the fourth tern, (ext ) Satiseto 
indicate that this term changes its sign such that friction 
torque resists spindle motion. This term will take on the 


same sign as the viscous friction tern, Oe ws ; 


43. The Pitch Setting Mechanism 

Two designs for the pitch setting mechanism are investi- 
gated. Both are hydralic. Design I is a pump controlled 
system using a variable displacement pump and a control rod 
running through the propeller shaft and is shown schematically 
in Figure A-7. Design II is a valve controlled system with the 
power piston within the hub. Design II includes a safety spring 
to make the propeller fail safe in case of loss of hydralic 
pressure and is shown in Figure A-8. 

The relationship between pitch and control rod or piston 
position depends on the particular hub mechanism involved. 
Several designs are in use. Figures A-8 and A+-7 show two 
possible methods. Figure A-7 shows two additional methods 
also in use. 

For a given hub diameter the parameter R in each mechanism 


will be approximately the same. The ratio A. for a given pitch 
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FIGURE A-=-7 


Design I Schematic 


=) Sinner 


a lic Force Amp. 





et, Spiny | (Ot. ‘. fy |e —~9 
Pilot Rod 
G \opucharse YS aa 
Lis Aiea Or 


LLL 


FIGURE A-8 
Design II Schematic 
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ratio ”/D will be only slightly different for any of the 
four mechanisms for small values of ie The greatest 
difference will occur for large /D. ‘Table A-1 shows 
the value of Kc/R for a ?/D = 1.5, which 1ls as large a 


pitch ratio usually encountered. 


Design I Design II Fig. 9-A(a) Fig 9-A(b) 
478 430 466 ES) 
Table Awl 


Values of the ratio %2 for a pitch ratio ?/p= Leo 


Since the mechanism of Figure 9-A(a) yields the simplest 
relationship between pitch and control rod displacement, it 
is used throughout the study. The overall results, however, 
will apply with reasonable accuracy to any of the hub 


mechanisms shown. 


ae (39) 


D R 


For Design I, 


K 
Kpx = fag +P, Sk +f, 2 (40) 
where 
x = pump control displacement, in. 


Kp = pump constant, jn?/in 


Kpx = flow from the pump, in? 


X, = control rod displacement, in 
Xe = velocity of control rod, in/sec 
Ap = control rod piston area, ine 


-4/7/- 








FIGURE A-9 
Two Additional Hub Mechanisms 
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FIGURE A-10 


Design I Block Diagram 
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XoAp = piston flow, in?/sec 





Po" = control red foree, 1b. 
K; = leakage constant of pump, in?/16 
F, Sle leakage flow 
AR 
Fe = time derivative of control rod force, 1b/sec 
Ke = v = compliance constant of the system under 
pressure, 1in°/lb 
= volume of o11 under pressure, in? 
4 

B = = bulk modulus of hydralic fluid, in 


2 2n 
= hydralic pressure, 1b/in@ 


oy 


Ie 
P = compliance flow, in?/sec ‘ 


Ke 
AR 


Equation (40) neglects leakage around the control rod 
piston and expansion of the hydralic piping. These factors 
can easily be accounted for increasing Ky to account for 
piston leakage, and increasing Ko to account for piping 
expansion. The compression of entrained air bubbles in the 
hydralic fluid can also be accounted for by increasing K, 
appropriately. 

For Design I, 

Fog = & Qz (41) 


where Fog is that component of control rod force transmitted 
from spindle torque, lb. 


Using equations (39), (40), and (41) the expression 
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for the dynamic motion of the pump-propeller system is 


Input Ram Flow : Leakage Flow x 
Kp*% = pesAe + a ae bas 





fe TL AN 
(fet ES) + (ah eR) 


+ he (BaheM.b) - See. ae?) * a ol (2g 


rae low 
re $rTN 
Ean Rye ae 


+ & (“aN (“ah tApR) + (Zea + MER) § 


(42) 

where R = crank radius, ft. 

fp = control rod ram sliding friction, lb 

Np = control rod and piston, viscous friction 
coefficient, Lb-sec 
ft 
Mc = control rod effective mass, 1lb-secé 
ft 
For Design II, as shown in Figure A-8 the equation for 


flow to the power piston is: 


(yexo) L Ky - XQAp + PAL + 1 ec : (43) 
AR AR 


Where 

y = pilot piston displacement, in 
L = inlet port width 

K, = valve constant, in?/in 


The control rod force in this case is the same as with 


Design I (the values of parameters will vary) with the exception 


=u) - 





of an additional term contributed by the safety spring. 
The resulting equation for Design II is: 





Mechanical 
Input Ram Flow Feedback Leakage Flow 
= / A Bxiulliz (« Ku) 2 
yl Ky AcOR mn Kx ihe es sow 
a Kev W m (Ko+ Kio) me, f+ fy a 
ie dR 
fn 
wl Sgt ApR) + 2 (5 Tia + MR) 
ahs reym Cow alPine ce Clow 
1 A RLe- -£) | + 1 A. {Cea En) nut 
ee 
_ KeVow ( Se par ear 
— Ko : + a r way) 
A Spring term 


KeQN | 
a 5 (G04 8) PH 


(44 ) 


Where 


> 


spring constant, lb/in 


spring force when 2 =O, lip 


ape 
ll 


The pilot valve-stroke piston system shown in Figure 
A-/ for Design I is standard and is fully analyzed in 
Reference 33. It can easily be shown that the time constant 
of such a system, loaded with 0.4 1b-sec* mass and supplied 
with 1000 1b/in® pressure oil will have time constants less 


than 0.1 sec even with excessive volumes of oil under 
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pressure. Since the transients of this portion of the 
system decay in a small fraction of the remainder of the 
system transients, the pilot-valve-stroke piston system 
Can be viewed as a simple batt. 

The feedback provided by the mechanical linkages of 
Design I are shown in Figure A-10. 

From Figure A-10 it is evident that 


maT, “He. | (45) 


Substituting equation (45) into (42) $ using relation 
(39), 


Input Ram Flow Mechanical Feedback 
KG, te oe _ 2Re Ke Ge 
re : ” 
C+, a (e+4,)IrD 





cakage ~ 
+ Kv a el Kav, w 7 (Kot Kio) . 
R R ic 


* (L.+ st), A ( fN 1h, R) 
i. & ( Zan a MeR) | Beam Na 


: ay: ae K,,] 
ale z= sa + (Cp ma | 
r é. (“se +Anr)t 





(46) 
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With the exception of the extra terms envolving 
the spring force in equation (44) equations (46) and 
(44) are of exactly the same form. They are arranged 
so that the analogue can be easily drawn. The analogue 
diagram for Design I is obtained from the analogue 
diagram of Design II by simply removing the two spring 
force feedback loops and adjusting parameters to match 


the actual systen. 
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APPENDIX B 


Determination of System Parameters 


l. Ship Characteristics 

The Victory Ship (AP2) was chosen as a prototype to 
obtain representative parameters for the system. The Troost 
B4-55 propeller curves were used to determine the propeller 
operating point for the Victory Ship specifications. Values 
for wake fraction were taken from references 16 and 17. 

These values are summarized as follows; 

Sip, vietory, AP2 

SEME | cues oc <0 OTe OO a cibls ec ccs so cee 6s OOOW pre 

Deepropel ler Glameeer . ..cccccccsececslOsed ft. 

STG ste |) ES 6 ccs ccc ccs cescenecss cook.) Kkmits 

eG, “He RAMOOl...« © aie eee 4 00's 00 6 0-00 ere « oo LOOneapm 

Number of propeller blades, N.........4 

Effective horsepower, EHP ....e.eeeee--4050 hp. 

Tee WOM eerie, «6 sass ve s0% sees OSOROOO Ube 

Propeller Torque, Qp..--sececeneceeeesdeLXlO? lb-ft 

THrUst. Walve ei ractiOn we ses ss eee 40 

Torque wake fraction,ug MAG e « « clés + RO DOO 

Thrust cOleifici@mt, K+ ..cscseccccceeOeoLO 

Torque coefficient, Kg pig 6:4 0 ove! ow 9 oo OE 


Pitemeretio, F/D .ocsscocsccccscowcces lms 
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Mean width ratio, MWR ...... 6 GbieeS oe oe Oe 


Speed Of MAVENGe, VO ..cccccccsccccsceelLO.! ft/ase 


2. Propeller Characteristics 

The Troost curves for Kg were then linearized about the 
operating point. The resulting linear approximation to these 
curves are shown in Figure Bel and Bee. 

Values for K,, Ky; Kes Ky» Kes and Ke were then obtained 
by substitution of the appropriate values into equations (7), 
(8), (9), (14), (15), and (16) respectively. 


- (1,024) (18.24) * _ 
= ATT, bg = 4O Tre lb-sec* = 288 lb-sec? (7) 


h 
K, AN eis — (.476) (1,024) (18.24 1pfege2 


2TT bg 20 IT 


863 lb-sec“ (8) 


5 
(AEE 024) (18. ° ) 
Hopp? 22 2 Lbesec’-ft 


Kz = 
: 4m bo 40TT © 
= 1840 lb-sec“=ft (9) 
3 2 2 
Ky = pp f (1...024) (18.25) ipoaee 
= 69.5 lb-sec@-ft (14) 


—-55 _ 


saaany OY pazoaUln 


Ql 6° 8° 7 dg Q° Se 


SE'O + “Nol +L(92b) = = 
y+ Dy B+ chu = F 
| 


| 
| 


Ge'0 =% 





b/ 


ay 





SOAANTD +y) Pezl4vaul4 


Ga sao io 


ae 3 = 
aa 6 Aa Z “a a i e a je ae 
Oo 
l'9 
t'9 
“3 

g £'0 
bo 
Bor AS c+ oO ~ 7 

- S 
"'a-= "4 ; i 0 
Lytdyty + piu = Ge 7 

ee 
“Y= OE'Z = =a aa ge, GO 
d VY 392° 60 
Oe, 

14"? 01 
i) 
| 
| Ay 
SLO'O = Ww aa Sy 
$'| 


Va 





i = Me D? _ a (.975) (1.024) (18.25)? Vib wee 


2T by 27 (2.5) - 


420 Lb-sec® (15) 


ft 


4 4 
Ke = hrpD* _ esr eee 1 pe amice 


ATT" b_ 


125.2 lb-sec@ (16) 


Curves of Kg were not available for any propeller. 
Rupp's data for control rod forces on a controllable pitch 
propeller while operating on a tug (YTB-502) in free route 
were analyzed as follows:* 

An approximate wake fraction was computed from a formula 
relating the block coefficient of the tug to the wake 
fraction. The propeller speed of advance, Vo, was then 


computed for the range of speeds tested. The speed of 


4, Lewis A. Rupp, "Controllable Pitch Propellers," 
Transations of the Society of Naval Architects and 
Marine Engineers, Vol. 56, 1948. pp. 272. 


5. S. Curtis Powell and W.V. Bassett, "Practical Aspects 
of Torsional Vibration in Marine Geared-Turbine 
Propulsion Units," Society of Naval Architects and 
Marine Engineers, March 1944. 
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advance, J, was then computed from its definition and the 
available data for each pitch ratio and speed. The correg- 
ponding control rod force was taken from curves given in 
Rupp's paper. Since these forces represent both the 
hydrodynamic and aentrifugal components of torque converted 
to force, the centrifugal force component of torque was 
subtracted from the total torque obtained from this force. 
The blade torque coefficient for each pitch ratio and 
advance coefficient was then computed from its definition, 
equation (19). The results were then plotted and found to 
be very nearly linear as shown in Figure B-3. The tabula- 
tion of the calculations leading to the curves of Ke are 
shown in Table Bel. 

As stated in Appendix A, the blade hydrodynamic spindle 
torque is a strong function of spindle position. The actual 
values for Kg could be shifted up and down the family of Kg 
curves by shifting this spindle axis. For the purpose of 
this study the Kg family of curves were assumed to be 
oriented such that the curve for zero torque intersected 
the origin of the pitch ratio - advance coefficient plot. 
This can be justified only in that a propeller could be 
made to give this result. This design would give a much 
more reasonable torque in the region of normal operation 
than merely using the tug propeller curves. The tug 
propeller was designed for a pitch ratio of about 0.7 in 
lieu of the Victory Ship with a P/D = it5. Whe eilore 
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FIGURE B-3 


Derived Kp Curves for YTB-502 
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Tabulation of YTBSO2 Free Route Trial 


TABLE Bel 


Calculations Leading to Kp Curves 


- 332° x 55 


(90.64) (25)(9.72) _ 


Voft/sec p=5.8 p=0.7 p=7.6 


1.05 
8.46 
9.89 
11.30 
12.0 
14.10 
Leo 
16.92 
Main BS 


82 ~ < 
100 89 80 
Lue 103 92 
136 120 108 
155 ILENE 23 
76 D2 140 
200 15 160 

- 200 182 

y= Yo 
nD 

544 ~ - 

529 ~=. 606 679 

524 86.595 661 

516 7905 652 
Os eect 635 
491 ~ 560 614 

- ° 536 588 

Je co 
S23, 
FIGURE B-5 


YTB-502 Hub Mechanism 
-6\|- 


~ 0.527 = Bloek coefficient 


Ca 0.10 = 0.164 = wake fraction 
2 





F, x 107 lbs (includes 
centrifugal force component) 


DO ell wel eel oo 
rOnNMrFWLO! 
DOOLAYHWOO 


800 


Ur NOVO FO 


al ool all ae 
PUFF OWUWL ! 


89 


a aell orll a 
NJUHOWOW! ! 


Qpy ft-lb-x107", centrifugal 
force component subtracted 


0238 


~O1L922 
eO1550 
01340 
sl 207 
01083 
~00910 


ANE EL 
e 6° e 
OVUINM FOW! & 


OM COU UI 





~300 


0302 
~0253 


~01740 
01354 
~O1112 
00930 
00820 


Note: / =1.024 


to match 
Troost's 
curves 





(m,,) and spacing (bp) were assumed to remain the same. 
The resulting linearized approximation is shown in Figure 
B=4, 


Values for Ka; Kg, and K. were then obtained by 


9 
substitution of the appropriate constants into equations 
(22), (23), and (24) as follows: 


m 4 2 
So Ce Qep te eo) lb-sec 


4TT?b. ft 


4 4 
ee m D - (.650) (1.024) (18.25) "_ 1959 lb-sec (23) 


2 Top 211(6) 


K 


HBPD? _ (0)(1.024)(18.25)? _ o (24) 
41 “bp ATT? (6) 


Sliding friction, as represented by equation (25), 
contains the thrust perpendicular to the spindle axis. This 
is the resultant of force vectors in the plane of the 
propeller disc as well as perpendicular to it. These 
vectors change with pitch, speed of advance, and RPM. These 
dependencies are difficult to account for, and Ty, was 
assumed equal to ‘To and constant. The dimensions assumed 
for bearing length and diameter are based on dimensions 
scaled from existing CCP hub mechanisms as shown in Ref. ll. 
In each case the dimension was assumed to be a value giving 


the least final value for sliding friction force. A friction 





SeAany (1 oUhposphy) Sy pownssy 
eel etal e' 


St 


- ' Ss’ 





uit isodl pesyo 
Of, Sepeq 240404 
OL Spuat e% 


e/ 


S‘/ 


Ya 
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coefficient, f£, of 0.10, was assumed. 


4 4 
Quo = fr LAN _ (0.10) (1.14) (9.8x10'") (4.56) (25) 
b (1.4) 


1.455x10° lb-ft. 


QBs 

Since Qpg is so large, even with low estimates for f, 
r, l,and |, and a high estimate b, the value for sliding 
friction contributed by the control rod and piston was 
assumed zero. i.e. fp = 0 

A sketch of the hub mechanism with dimensions is shown 
in Figure B-6, 

Reference licontains a plot of control rod force for 
the transient conditions of going from full ahead to full 
astern while the propeller was spinning in air. This curve 
is shown in Figure B-7. For a spin test in air the thrust 
forces causing high values of sliding friction are small. 
The long region of both transient curves that are parallel 
with the steady state (centrifugal component only) corres- 
pond to the condition of maximum pump output and constant 
speed of the piston. The deviation from the steady state 
curve in both cases is approximately 3 x 10* lbs. 

Reference ll states that the total control rod travel is 
0.318 ft. and that complete reversal of pitch takes five 
seconds when the propeller is spinning in air. 

Assuming that this entire force is caused by viscous 


friction, and the maximum ram speed is equal to the 
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FIGURE B-6 


Propeller Mechanisms Showing Only One Blade 
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& = Control rod force in lbs - Thousonals 


45 


Emergency ahead trom 


full astern 


40 

3 

30 

2 

S 30,000 /b, 

1S 

10 enter p 

5 £29] = sy , 

0 S¥a7. 







30,000 /6. 


Emergency 
astern from full ahead 


oS gy 6€u5 6% 3 21 0 t 239 4+ 576 7 eee 


<—Ahead Pitch in Feet 


FIGURE B-7 


YTB-502 Control Rod Force for Spin Test 
in Air from Data by Rupp 
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average speed, 


X. = on = .0636 ft/sec 


the viscous force = 5.0 x 10#L> = 4/ z. 


i. 5 lb-sec 
A = 4,71 x 10 ae 


Boer ey — total viscous damping coefficient, oe 


Assuming the viscous forces scale as the area and 


using the viscous friction term of equation (42), 


1 (Asn AR = Ce oR (47) 
TD oY (2 YTB-50° 


i N 
— B _ 4 
WD = + ARR) = 4.71 x 10° (28.25 )? 1,2 


CaN + ARk) = 3.65 x 102 iLbssec 
R ft 


In order to obtain a value for the polar mass 
moment of inertia of the blade about its spindle axis, the 
blade is assumed to be an elliptical plate with an average 
thickness = : fe 


The formula for the moment of inertia of an ellipse 


about its major axis: 


BY 2 


where t = thickness of plate, ft. 


(48 ) 
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(> = density of metal, Lbssecr 

a = length of minor axis, ft. 

b = length of major axis, ft. 

Japp = mass moment of inertia of blade about major 


axis, lb-ft-sec 


The axis of the spindle is offset so that approxi- 
mately 40% of the area is forward of the spindle axis. 


3/2 
° — a bil N e+ 
q J4 pu = tP, =— + © me abN 


where 
C = distance of spindle axis from the major axis, ft. 
N = number of blades 


Substituting values: 


y (pope, (5e57) (@s125) 4 
B Soe) “aaa (49) 


5 Gp 2h (eee ) TT (5.37)(9.125)4 
6 32.2 


Jy pil 


808.5 lb-ft-sec* 


Jqph 


Estimating the mass of the control rod a length of 20 


feet was assumed. 


Mo = PoTl (Re )2 1 = $22 (2.65 )2 (20) = 740 Lb=sec® 
12 Sane 12 ft 


Po = mass density of steel, 1b-sec® 


ft 
radius of control rod, 7,67 inches 


9) 
iH 
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L = rod length, ft. 

This is a large estimate for control rod mass, but it 
is shown in Appendix C that inertia loading can be neglected 
even for this large estimate. 

Allowing a 30% increase for entrained water an estimate 
of J4p = 1050.0 lb-ftsec* is obtained. 


The assumption of an elliptical plate was also made in 


evaluating the constants for centrifugal blade torque, Kio 


and Kj. 
The formula for an ellipse: 
2 2 
x4 =) 
aa Ra . (50) 


Translating the axis to the blade root and solving for X: 


2 
rv +2 (51) 


(~ - 
The blade section area, A, is a function of distance from 


the root, 


2 
betormar | aG- Ete 2 | ; (52) 


The distance from the center of the area to the spindle 
axis in the x direction is assumed constant = c. The 
distance of the center of the blade section to the spindle 
axis in the y direction is a function of rake. 

y = n,/f (53) 


where m, = slope of blade from vertical. 
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Substituting equations (52) and (53) into equation (36): 


tip 

Kio = NC, \ axyaf 
ore: 

K,,=Nnf b 

10 P (2 at | ata - 2 - + eB om. aff 
hub 


If m, is small the vertical distance is approximately 
equal to the blade length. 


tip 
Kip = N 2 \= E (4+ = r 3 + B | cm,YdaY . 


hub 


Integrating: 


ee 

NW 

ol 
Lae | 
re) Milo’ 


+ N Pp4tem, 5 z 
2 


NM 
Seal 
oO Nfo’ 


Assuming m, = ¢ and substituting values for the 


remaining variables: 


Ko = 4 (222) (S ) (.67)(4) (2.685) 24257 


Jee (2) (8) (4) 
(9.125) 
4 (aude +)(.67) (4) —— 
: wt % ; (64) 
ges hese Lb-ft-sece 


(36) 


(54) 





The formula for the moment of inertia of a rectangle 


about it centerline: 


ta BX 
MIN 12 : 
where the length = X (6S) 
also 
eet 
ee = Sy: (56) 


Using formula (51) for X, the expression for Kj, can 
be approximated using equation (37). 


re 
Ki, = BTID e, \ ux - Inry) oY. (Si) 


“ne 


b 
«lhl 2 [aca- 2 2 7 py)" 


: L a(d - + Seay} (57) 


After expanding the cubed term, (57) integrates to 


a 


_ NV D 
Ki, = NTT Df tab L T +2 + a2 | 
12 64 — 
32 
~NITDP ta BY | _ NWDe ¢3p2 
Fe “a ] 32 a (58) 
2 


Expression (58) reduces to 


= NT D 

Se NT DE, bt ae | 22+ (602-3) | -sent2} (59) 
768a 

Upon substitution of the limits. 


Substituting values for the various parameters: 


eo — 
Ki, = ATT (18.25) >) (9.125) (4) (5.397 6.37 
goo) (Soon) 


2 
+ Tr 6(9.125)° = sm | -32 9.225140 
Ky. = 945 Lb-sec* 


3. Pitch Control Mechanism Parameters. 

The value for R is obtained from the geometry of the 
hub mechanism and is shown in Figure B-6 to be about 1.2 ft., 
which is the maximum feasible for this size propeller. 

A maximum operating pressure of 1000 (b/in® was assumed. 

The piston area required is obtained from summing the 
forces on the control rod and dividing by the pressure. 

Combining equations (21) and (35) and substituting 
Operating point valves for the variables yields the spindle 


torque cause by hydrodynamic and centrifugal forces. 


QpH + Qgq = (K7 - Kyy) pWe - KgVye ~ (Kg + Kyo) © 
Qpp + Qpq = (482- 945)(21.0) (297 )e 

z 1955 (16.7)(201) ~ (0 + 105.7) (LOM ) 
QpH + Qpq = 1.063 x 10! 1b-ft 
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when the maximum torque condition prevails. 


ore f trt, 4 nN = (.10)(1.14)(9.8x10") (4.56)4 
D 1.4 
Qpg = = 1.455 x 10° id-ft 


If these were the only torques to overcome, ie 
neglecting viscous friction and inertia loading, the total 


torque would be the sum of the torques above. 


QBH + Qng + Qpg = 1.063 x 107 + 1.455 x 10° 
Qpy + QBc + Qpg = 12.09 x 10° 1b-ft 
Dividing by R yields control rod force. 


_ Spy * 8p¢ * Sas = 12.09 x 10° lb = 1.005 x 10! iy 
. R 2 


The ram area, Ap, becomes 


1.005 x 10! 1b-in® 


1000 1b 


ae Geer © By ¢ 10+ ine 


A piston of this size would not fit within the hub of 
4 ft. diameter for the Design II type mechanism and is 
excessively large in any case. Hence a piston area of 7/75 
in? is assumed which gives a piston radius of 15.75 in. This 
size is about the maximum that can be housed within the hub 
and for either design there is a likelihood of stalling 
under the heavier load conditions. <A larger diameter 


causes extreme demands on pump output in order to obtain a 
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reasonable response for the pitch setting mechanism. This 
pump must be very large even with this size piston. 
As is shown in Appendix E, the time constant of the 


pitch control mechanism is 


Tp = RAp 
Kp TT DM 


(60) 





It is necessary to reduce this time constant to as small 
a value as possible. The larger this time constant, the less 
effective the control mechanism becomes in RPM regulation. 
This is shown graphically in Figure E-6. As Tp increases the 
frequency at which the ratio of RPM fluctuation to wake 
disturbance begins to increase for all systems with regu- 
lation decreases. 

Higher values of Ky result in lower values éf 7. The 
higher values of Kp correspond to large pump capacities. 
Assuming a 300 H.P. variable displacement is the maximum 
available Kp is calculated as follows: 

A 300 HP variable displacement pump is approximately 
capable of a maximum delivery of 120,000 in?’/min = 2000 
in?/sec. 

The maximum control arm displacement is arbitrary and 
depends on the length of the lever. A reasonable value 


for this parameter is 4.368 inches. 


Hence, 
= 2000 in? = = 458 in? 
e 4,368 inesec in<sec 


BS 





The leakage constant, Ky,» can be estimated by 


assuming a maximum % of full pump displacement at full 





pressure. 
pee 
Ki = S—- (assumed %) 
L000 a 
in 
K. = 0.2 in? for 10% leakage 
L spec=-lb 
in? 
Kr = 0.02 for 1% leakage 
sece-lLb 


The mechanical feedback shown in Figures A-7 and A-10 
is determined by the length of linkages 6 and 4% or the 
length and width of the valve port in the case of Design 
II. Increasing this feedback decreases Tp, but also 
decreases the overall pitch mechanism feedback gain as 
shown in Figures E-4 and E=5. This corresponds to requir- 
ing a progressively longer travel of mechanical input at 
y in Figure A-7 or on A-8. The sensing devices whether 
tachometers or torquemeters will produce electrical or 
weak mechanical outputs. A long mechanical travel at y 
in either case introduces additional time lags that reduce 
system performance because of inertia loading at a point 
in the circuit havingalow power level. Rksebileabie values 
for the linkage lengths were therefore assumed as follows: 

e = 5 inches 


f = 15 inches 
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5 inches 


~ 
tt 


ox 
i 


5 inches 
The Go which corresponds to the gain of the stroke 
piston system is assumed to be one, i.e. it acts purely 
as a force amplifier. 
Wz Leese = dlejil lee = 0.063 
(e + f )1TD (5415) T (18.25) 
The time constant for the stroke piston system is 


smaller than any other time constants in the system by 


several orders of magnitude and is therefore neglected.° 
The compliance constant, Ky is determined from 
equation (62). 
V 
kK. = = (62) 
as 
where 
V = volume of oil under pressure, in? 
2 
B = bulk modulus of hydralic fluid, — 
B = 0.27 x 10° =, for standard hydralic oils. 
n 


The maximum piston travel is, 


Xo = Rp = (ame) (21) - 0.44 ft. 


Max WD TT (18,25) 


The volume of oil under compression in the piston is, 





aye Truxal, JeGes O Bean p.15-44,. 





AR Xo.) = Oot# x 775 x 12 = 8,200 in? 


Assuming a 10 foot length of piping under compression 
with a diameter of 3.4 inches, the volume under compression 
in piping is, 

(20) (12) (3.4) "ry = 1180 in? 


The total volume under compression: 
8,200 + 1,180 = 9,380 in? 


5 
, V _ 9380 = in- 
ek. == = = BOS 50 
B” 0,27 x 10° 16 


4. Power Plant Characteristics 

Propeller inertia (including entrained water), 
propeller damping, power plant inertia referred to the 
propeller shaft, and engine damping assuming the turbine 
is a constant torque device, is taken from Sutton's 


vibration analysis of the Victory Ship, (apo). ! 


3.09 x 107 lb-ft-sec* 


| 
Jip 


Jip =i oe 10° lb-ft-sec® 


ma 4 


ft-lb-sec 
rad 





7. A.D. Sutton, "Torsional Vibration of a Geared Turbine 
Propulsion System." New England Section Paper of the 
Society of Naval Architects and Marine Engineers, 
October 1953, p.27. 
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oa 7.54 x 10° ft-lb-sec 
iP ‘ —— 
rad 


* 
Ap is also found in Appendix E to be equal to Kz 


and can also be determined from propeller curves. 
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APPENDIX C 


Development of the Non-Linear Analogue 


The final expressions derived in Appendix A could 
be programed directly for analogue computation. Simpli- 
fication of equation (46) Paatiy reduces the complexity 
of the program and is therefore desirable. 


The compliance term in equation (46) contains 


oS (Ro-Eia) pw2 - X8Vow =~ (Kg +Kio) yy2 
at R R R 
t (ep + fr HX OM) (63) 
bR 


Performing the derivative, 


ih ‘alae y 2p Ww eae KgV - Kg WV 
R 
+ O. (fp + ert, XN ; 
dt bp 


i is the variable being regulated, and even without 
regulation W is small. Therefore, those terms containing 


c) are assumed zero. 


The term (fp + £7, L N_ ) is pictured as a function 
19) 


of p in Figure C-1l. 
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FIGURE C-l 


Sliding Friction as a Function of p 
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This sliding friction term is discontinuous at f = QO. 
The time derivative of the friction term is infinite as the 
direction of motion of pitch is reversed and therefore 
causes an impulse of flow or a shock. It does not allow 
a volume to flow and therefore can be neglected in the 
analogue, The presence of the term indicates that a 
hydralic shock will be generated for each pitch reversal 
in the actual system. 


Expression (63) reduces to 
{ A 
R 1 (K7 - Ky) pw® - Keto an 


with the above simplications. 


Expression (46) contains two terms in f. They are 





TID Ap j 
and 
Bo ey og ur ) Ke 
& a , 


Substituting values into these two terms: 


D en + TH0(2.2) 0350 - 3 (52.5) 


1D Wa AR TT DAp 

and 
Bb (2.08 x 107) 2, =b (4,16 x 10%) . 
TD “R DAp 





Hence the first term, envolving inertia and leakage 


can be neglected. 
Substituting expression (47) for expression (46) in 
equation (37) results in the following three terms 


containing p. 


- oR ~ Wise Meo) = — jou, 3 
Pp AR TT (18.25) af = 
8 Gh ee (‘Se +Agk ) = (,2)(1.742 x 10°) : 
TT DAR R (18.25) (775) 
= 0.785 p 
(For 10% leakage) 
p Kel&rRiy) ww? (.0336)(480-945) 2ST) 


Ap 775 


= 2,20 p in?/sec 


Hence, ignoring these last two terms, equation (46) 
reduces to: 


R Ky, |- 
KA = p TD AR : 7. | (K7-K,}) pw® ah KeYo W 
R R R 
-(Ko + Kio) aw faerie frt, 2 N | (65) 
R bR 


AR 


Ke 
R is Tei aaa 


= SV 





Equation (65) and the remaining equations for Design 
I derived in Appendix A are programed for analogue compu- 
tation as shown in Figure C-2. Substitution of the 
parameters determined in Appendix B is shown in Figure C-3. 

The analogue computer used will not accept signals of 
greater than 50 volts. Also, when the output of any 
amplifier exceeds 50 volts the amplifier response becomes 
non=linear and saturates completely at about 57 volts. 
For this reason it was necessary to use the following 
scales for variables: 


pitch, 1 ft = 0.1 vélts 
shaft speed, l rad _ 0.4 volts 
sec 


speed of advance, 13 = 0.1 velts 


A scale larger than this saturates the path containing 
pu terms. 

Care must be taken when using different scales for 
variables when non linear elements such as the multipliers 
are used. 

The general rule for multipliers is to divide the 
amplification following the multiplier by the product of 
the two scale factors of the input variables. 

Scaling through linear elements such as amplifiers 
and integrators is accomplished by increasing the coefficient 
for each input path containing pitch and speed of advance 


as a variable by a factor of four. 
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Finally it is necessary to move various gains around 
the system, following the rules of block diagram trans- 
formation, so that individual amplifiers will not 
saturate and the inputs to each amplifier do not exceed 
four, all with coefficients having their decimal point 
within three digits of each other. An example of this 
last step is shown in Figure 0-4. The gains must also 
be arranged so that no amplifier is called upon to 
provide a gain greater than 10°, 

The portion of the analogue marked "Selector- 
Amplifier Section" effectively produces the response 
pictured in Figure C~l, The bias and amplifier gains 
are set for a very small sinusoidal input to give the 
correct voltages corresponding to the term, (fp * £,tl w ) ; 

DR 
Figure C-6 is the analogue with all scaling 


completed and the program ready for the computer. 
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FIGURE C-4 


Example of a Block Diagram Transformation 
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APPENDIX D 
Non-Linear Analogue Oomputation Results 


1. Difficulties. 

The Philbrick Analogue Computer of the Mechanical 
Department Systems Project Laboratory at M.I.T. was used 
to obtain data. This computer utilizes electronic 
multipliers to perform multiplication and division. The 
multipliers produced a hum that was within the band pass 
of the amplifiers. Although correct average values for 
the variables, WW; Vo, Vea, and ear were Obtained, the 
hum superimposed on the signal was transmitted through 
each amplifier and in some cases obscured the data. The 
power plant-propeller system and the pitch control 
mechanism were operated separately with satisfactory results, 
but the connecting link, marked G in Figure C-2 introduced 
hum for the variable p saturating the integrator that 
followed. For this reason, no overall system response was 


obtained, 


2. Useful Results. 
The response of the power plant-propeller system to a 
sinusoidal input of V, was obtained. A two channel Sanborn 


recorder was used with the input, Vp, on one channel and the 
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& response on the other. The Sanborn recorder stylus 
deflection was calibrated against a computer amplifier. 
Each amplifier has a provision for adding up to * 50.0 
volts to the output in increments of 0.1 volts. With 

one volt per cm. deflection the recorder output was 

linear over * 2 cm. The maximum output recorded was 

+ 2.2 cm. and the accuracy of readings is approximately 
O.l1 cm. The results are shown in Table Del and are plotted 
in Figure E-6. 

The time constant of the power plant-propeller 
system was observed to be approximately 3.28 sec. This was 
determined by injecting a step of Q and observing the 
response in W, The decay to steady state in ™ occured 
in about 10 sec. or three time constants. 

The values of Qp and Qp generated had an average 
value within 5% (3.1 volts and 7.0 volts respectively) of 
the calculated values at the operating point, but a sinu- 
soidal hum of volts was superimposed on both, 

The node corresponding to the variable pf had a sinu- 
soidal amplitude of .6 volts for a 8 volt sinusoidal input 
for Vg with no friction feedback. With friction added this 
signal was reduced to less 0.1 volts (the approximate level 
of noise). This is an indication of the importance of the 
strong stabilizing effect of sliding friction. 

The feedforward path containing 7, had no effect on 


p even with an input of 30 volts for AVo. 
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TABLE D-1l 


Power Plant-Propeller System Frequency Response 


frequency AVo 
cycles/sec Input wie, 1Olog |Aw 
(Volts) (Volt Sx¢) |A = LOAV , 
O B08, Cals 0.050 ae L329 
eO1L 3,0 O.15 0.050 - 13.0 
02 3.0 0.15 0.050 - 13.0 
03 3.0 0.20 0.066 - 21.6 
04 5.9 0.20 0.066 = P18 
-(0}5) 3.0 0.20 0.066 - 11.8 
06 550 0.20 0.066 - 11.8 
~O7 Sio O220 0.066 = je 
08 570 0.20 0.066 = TiS 
09 4.0 0.20 0.050 - 13.0 
ow 4.0 0.20 OFO50 = ED ae 
2 8.0 0.38 0.048 - 13.2 
23 LiGO 0.62 0.056 - 13.4 
4 15.0 0,84 0.056 - 13.4 
35) 1e.0 0.79 0,044 = (8355 
M6) 24.0 1g Ol. 0.043 - 13.6 
of 26.0 Polo 0.042 = oon 
=0 28.0 Weis 0.0415 - 13.8 
J 50.0 1.25 0.0415 - 13.8 
0 BLS Ae, oe 0.0400 - 14.0 
2.0 5830 1.45 0.0380 - 14,2 
5:50 38.0 1.0 0.0264 Ee =, 
4.0 38.0 O50 0.0130 - 18.8 
5.0 40.0 0.00 0 - OO 
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APPENDIX E 


Linear Incremental Analysis 


1. Simplification of the model equation (46) contains the 


tern, 


12Gn e 
Sears @, 


(e + fy \rD 


which is mechanical feedback. 
Define 12G5eR = (66) 
(e + fo )7D 


As shown in Appendix C, a reasonable value for this 
term is 


M = 0.063 ; 


{TD 


As shown in Appendix D, the terms P_ (CBN + AAR) Ke 
R AR 


Kr, Kev 
and oe oS ee » are small in comparison to the remaining 
Ap R 


terms of equation (45). Neglecting these terms the 
analogue for the pitch setting mechanism, as shown in 
Appendix B, can be reduced to that shown in Figure E-l. 
This is equivalent to neglecting the corresponding term of 


equation (46), which reduces to 





2 
Dp = _— ‘A Kpx - = J (K7-Kj1) — a = 
TD * “ 
~(Kg + Kio) ws (fp + frt, JN ) (67) 
R bR 
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FIGURE E-l 
Simplified Pitch Control Analogue 
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FIGURE E=- 


Pitch Mechamism Amalogue Neglecting Leakage 
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if kL is sufficiently small the leakage feedback 
AR 


shown in Figure E-1 can be neglected and the system reduces 
to that shown in Figure E-2,. 


The transfer function for the system shown in Figure 








E-2 is 
, KD 
S/ = 
yis Rip 
Ky TDM 
S + <P (68) 
RAR 
This oonstitutes a type zero system having a simple 
lag with time constant. © 
Tp =k 
Kp WW DM (69) 


Substituting values determined in Appendix C, 


T> = (1.2)(775) = 1563 sec. 
(458) 17 (18.25) (.063) 


Leakage has the effect of reducing this time constant. 
The two largest contributors are friction and the pitch 
component of blade torque. Figure E-3 shows these 
additional feedback paths added to the system shown in 
Figure E-2. 
2 


The system transfer function that includes the pw 


feedback term is 





8. J.G. Truxal, Control Engineers Handbook. The Maple 
Press Con 5 Yorks, Pa. 195 9 D.c=55. 





FIGURE E-3 


Pitch Mechanism Analogue With Major Leakage 
Factors Included 





KpTt D 





Es) = gh we 
wa 5 + K pM TT DRA 4 (K7=K ) TD w K ’ (70) 
ancy. Ag R 
Ap°Ro 
R 
Pp ~ yf Si ee a ee (71) 
T= (775)* (1.2)? 
P (458) (063) TT (18.25) (1.2) (775) + (276-482) a0 Wa25Ke) 
Tp = 0.541 sec. 


The computation above is based on a leakage of 10% of 
pump output at maximum flow. The leakage caused by friction 
can be approximated as reducing the pump output or effectively 
reducing Kp as shown diagramatically in Figure E-4. 


The effective Kp becomes: 


ct = ¥P™nax - oo) Ky 
7 a =R (72) 
XMAX 


1 _ (458)(4,368) = (1.210 x 10°) 38 
ss (4.368) 


Ky = 384 in? 
in-sec 


The new Tp becomes 


Pxbirpm (384) 1 (18.25) (0.63) 
= = 0.654 


— 9 /— 





Pump Out put “Kee 





Xm mM inches 


FIGURE E-4 
Effect of Leakage on K 





FIGURE E=5 


Regulation by ‘' Semsing 
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The bad effect of increasing Abs is shown in the gain 
plot of Figure E-6. As Tp increases, the gain defined as 
the ratio of shaft speed variation to make fluctuation 
increases and at progressively lower frequencies. It is 
expected that 10% leakage is excessive and the T, computed 
for 1% will be used in all future calculations. That is 

T, 


e 
Viewing propeller torque as a linear function of p, 


= 0.560 6 


W, and Vo, 
Qp = KiVo + K5p + Kew + Qro 
where 


Qpp = Operating point torque, ft-lb 





2 Q 
ky = — 
Vo PD »W 
a >Q 
Ky = aa: 
> Pp Vo. 
dQp 
3 P 





Values for the above K*'s are obtained from propeller 


curves by allowing the variable in the denominator of the 
partial derivatives vary over + .1 of their operating 
value, obtaining the Kg at each of the two points and 
solving for the Qp and dividing the result by 0.2. 


Calculations for Sif K5 and K5 using Troost's 


-9Y- 


B4-55 propeller follow: 


hd, 25 
lr Ad eae (79) 
AVo  (.2)nD 


AK, = Ky (for J = 1.lJ,) - Kg (for J = 0.93) 


where I is the operating point advance coefficient. 


10,9 5 
ae ie © (-.007) (1.024) (FF (18.25)? _ _¢erg Lb-sec , 


(12) (=e) (18.25) 
oP (.2) px D 
D 


where AKg = Ky (for f= 1.1 (B) ) 


& if P=o0O, p 
— D ‘ (5 ),) 


ee (,02) (1.024) (F° F(18.25)* = 2.96 x 104 1b. 
; (.2) (1.15) 


In this case p/D and V, remain constant 


= 
but J = Ua, changes as n = ce changes by - O.1. 
nD Pail 


ae = (006) (1.024) (A2)* (18.25)? = 7.54 x 10* ft-lb-sec 


3 (.2) (22) 2 TT rad. 
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In summary: 


-6610 lb-sec 


A 
i 


Ke = 2.96 x 10+ 1b 


3 rad. 


A 
i 


Using primes (‘) to indicate incremental quantities, 
equation (74) is incrementally, 


cy 


f * 
oF + Kz WwW (81) 


Rp -_ Kivi + K 
From equation (10), 


§ 


or using the Laplace transforn, 


a a] (C+ CL) + Ign |= KS p(s) + K¥V'(s) (83) 
ow (s) =_ K5P(s) + Ki Vo (s) (84) 


ee ' 
(Cn tC ole ins 
The block diagram for which is shown in Figure E-3 


ec. Shaft Speed Sensing 

Sensing angular velocity and adding the pitch setting 
mechanism to the diagram of Figure E-3 is shown in Figure 
E-5. 

Define: C,=L, +%, = 3.08x104+7.54x10%=10,62x10% Ltab-see 


rad. 


—/O| = 





v=. 91t = 3.48 x 10° 


ae —— = 5,20 Gee. 
~t on 10.62 x 10 


G = gain for angular speed sensing 


The transfer function between Ve. (8) and u) (s) is; 


w'(s) _ Ki (1- T,S) 
of 


Yo (8) (85) 
(T $41) (Tp+1) +%o Gus 


7M 

The selection of G,, is usually based on criteria for 
transient response or stability. As G,, increases the 
effectiveness of regulation improves for lower frequencies 
but more severe peaking occurs near resonance. A compromise 
value of G,, giving a damping factor, P= 0.4 was chosen and 
used as the criteria for and RQ. sensing designs. A Gy) 
giving a damping factor of one provides essentially no 
improvement of regulation overthe fixed pitch systen. 

The damping factor is defined for the transfer function 


when the denominator is placed in the form: 


S° + 2— Ws + We (86) 


where 
BS damping factor 


arn = natural undamped frequency, sec7l 


Placing the denominator of equation (85) in the form 


of expression (86), and solving for Gu) and Ww with f= O.4 


(OG 





yields: 


ag (Tet Tp) gM _ (3.284.560)? (3.08x104) (.063) 
— p = 


Te Tp( 64) KS (3.28) (.560)(.64) (2.962104) 


- 0,82 ft=sec (87) 
rad 


Dua EG wb ~ 2n06/re1o" (.82) 
cle Vp M (3-08 x 10%) (. 560) (3.28) (.063) 


= 2.61 seo (88 ) 


System response for angular speed sensing is plotted 


in Figure E-6,. 


3. Torque Sensing 

If a torquemeter were used in lieu of a tachometer for 
sensing, regulation of shaft torque might be accomplished. 
This in turn would regulate shaft angular speed. 

Figure E-7 shows the block diagram for the system using 
torque sensing. 


L 
The Vie to W transfer function is now 


Ke 
wis) _ ke (T7841) 
V6 SUE IPE UUUEEIEE “Sy aUEE EERE UREN GEER SREP 
s*+ apa ae P XoGgT s+/l4+e3 + K3GQ 
C+ —M TT M 
T Tp 


(89) 


Using the same criteria of f= O.4 in the setting of 


Gq leads to a quadratic equation in GQ. 


sy Eke 


1} 
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FIGURE E-7 


Block Diagram for Regulation by Torque Sensing 
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Solution for G, follows: 


Q 

From expression (89), 

Tp+T+K3 Tp + K3G9T 
— aL 


ee i =2u,F (90) 


.560 + 3.28 + + 1.54x10°(.56 Sa toe gon10"sg (3. 28) 


10,62x10 
(.563) (3.28) 


= 2(.4)w,, 


2.301 + 8.35 x 10°G9= BW, 


2.88 + 1.042 x 10°Gg =w,, (91) 
also, , 
aoe ae K5 + 36g im \ 4 (92) 
Lm T Tp 
w= 1 + 7.54 x 10% + 2.96 x 10%G9 8/1 | 2 (93) 
10.62 x 104 063 3.26x.560 
w = (1.708 + 4.70 x 10949 )# (94) 


Combining expressions (91) and (94), 


1.708 + 4.70 x 10°¢9= (2.88 + 1.042 x 10°G, )2 


Rearranging terms 


+ 1.09 x 10774 5.53 x 10°Gg+ 6.09 = 0 (95) 
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The positive root of this equation is approximately 
+ 6.78X\9! ft =Gg This value of Gg is not un- 
lb-ft 


reasonably small when it is noted that the steady state 


shaft torque = 3.1 x 10° lo-ft. 


Substituting the relevant positive value into equation 


(91) yields, 
6 


n 


UO = 2.17 seo7- 


Dye aioeC sore slo” 6.76. a 107! \usees! 


6 
the “ (8) is plotted in Figure E-6 for these values 


Vo(s) 


4, Thrust Sensing: 


Propeller thrust can also be considered as a linear 


function of V>,W, and p just as Qp. 
T= Kon + Ke p + K6 


where 


K, = or 


2 | at the operating point 


% OT 
Sp jat the operating point 


‘? 
1e2 
i 


+ DT 
S25 | at the operating point 


mA 
OV 
il 


(96) 


(97) 


(98) 


(99) 


3 * 
Values for Ky, , Ke » and Ke” are found exactly as were 


=fay 


values for previous stored variables. Again using Troost's 


Curves: 


10,2 
K,* ee. Kp 0p°D* I ~.04) (1.024) ‘e=) (18.25)* 


Av, _ (-2)nD (.2) (22) (18.25) 


K,* = -1.325 x 107 Lb-see , 


ft 
—— at = _Knbn2@p* = (08) (1,024) (49)* (18.25)" 
; is (2) D (oP isabella 2755), 
K, = 6,00x 107 1b 
a an 


and 


a 
D> 
a | 


mie AZ. KppneD* - (04) (1.024) (28)° (18,25)4 
: (.2)(n)(2T) (2) A8)277 





i 
[> 
= 


+ ay Ny, x 10% lb-sec 
rad 


mv 
OV 
iH) 


Incrementally, equation (96) becomes 


1 # ‘ 
T = Ky Vii+Ko p + Kg" 


(100 ) 


The block diagram for the system sensing thrust is shown 
in Figure E-8. 

In order to determine the effect on ™' of thrust 
sensing regulation it is necessary to add the paths that 
generated Uw from Nig D> and sag This is shown in Figure 


E-9. 
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FIGURE E-8 


System With Thrust Sensing 
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FIGURE E-9 


Block Diagram of Complete Thrust Sensing Regulation 
Showing Generation of 
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The transfer function between Vo(s) and w (s) is 


more complicated in this case due to the presence of two 


' - ik K4GpKo* | 
W_(s) - 3 Ky" + 


(T>S+1) 


paths and two loops. 


(101) 
+* 
Vo(s) clea Ks Gp m Ke*x *Gm 
M(VpS41) MK (TpS41) (TS41) 
Rearranging: 
Ay -( Lg Ko"Ky*e 
ws) = LTT p Le 1, (102) 
z.¢ | a 
Vo(s) ge (2 + 27 n KeG ayn the Gn 
TpT MTp fr NT p 
tt + 
= K6 Ko Gp 
MELT Tp 
Substituting Values: J thio 2.46 x10%7. 1325x1047 
66 0 paiiidadeneLieise ss 
ws) . Pers COD a TT (10,62) (7.28) (156)(.063)10* 
er ne a a a. eT aS 
VO) S5 + 2(3.28+.56 6x1074r a i 6 ¥t0 Ar 
ee Nar VEY 
(3.2 8) 56) 7,063 aa (8 (3.28Y§ ” (.06 3X. 56 (3128) 


_ 2 479x104 (2 abr )aT 
TO 
~ 1, 063 Xior br i0*K2 128) (:56) 


sO = 





§ 
Wy — = .01895 S = (-0339 - 3190 Gm) (103) 
To Se ES ee Fe sae arr 
Vols S +(2,09+1.7x102Gp)S + (0.545-7.6x10°Gm) 


Again setting § = 0.4 


2 fl, = 2.09 + 1.7 x 10? Ge 


O) , = 2.61 + 2.12 x 10? Gp 


Solving for Gp: 
(2.61 + 2.12 x 10°G,)? = (.545 - 7.6 x 10°Gn) 


10 


4.5 x 10° Gn* + 4.31 x 109G, + 5.26 = 0 (104) 


Equation (104) has only imaginary roots. In fact 
setting Gp=O in equation (103) results in 

Wo, = 0.736 sec7L 
and é = 1.42 , 


as well as a gain at zero frequency 


] 
@ (0) ~ ~ 0.0621 ; 10 log, |%_ (9) | = -12.06 
veo) v' (0) 








which checks with the fixed pitch case as it must. 

For all values of Gp greater than zero, é increases 
and zero frequency gain also increases. 

For negative values of Gm, the pitch control mechan- 
ism would sense a decrease in thrust and respond by decreasing 
pitch which would further decrease thrust and the system 
would settle at zero pitch, zero thrust, and maximum ™~ . 


Assuming a value of Gp = + 107? ft , the unfortunate 
Lb 


dilatorius effect of the thrust sensing can be shown. 


Sl - 





Equation (103) becomes 


ae = +.018955 + .0658 (105) 
Va\s s* + 3.79 S +(.469) 


Equation (105) is plotted in Figure E-6 for 


comparison purposes. 
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